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PURPOSE:  The  purpose  of  this  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  is 
to  document  techniques  for  computing  Factors  of  Safety  (FoS)  for  seepage-related  soil  instability 
of  earthen  structures,  using  the  simulation  results  from  a  3D,  finite-element  groundwater  model. 

BACKGROUND:  The  FoS  analysis  has  been  used  widely  in  geotechnical  practice  to  describe  a 
soil’s  state  versus  some  critical  condition,  often  associated  with  engineering  failure.  It  is  a  simple 
and  reliable  supplement  to  conventional  analysis  that  accommodates  the  combined  effects  of  the 
parameters  and  uncertainties  involved  in  the  calculations.  Traditional  practice  for  seepage  analysis 
has  been  to  use  2D  groundwater  models  (e.g.,  SEEP/W  (http://www. seo- slope. com/products/ 
seepw.asyx)  and  SEEP2D  (http://www.aauaveo.com/sofiware/2ms-seeD2d)),  to  evaluate  issues 
related  to  piping  and  seepage  through  impoundment  structures  such  as  levees,  dikes,  and  dams 
(Davis  et  al.  2009).  Although  the  incorporation  of  variability  in  model  parameters  into  these  FoS 
analyses  helps  to  account  for  uncertainties  from  different  sources  (Duncan  2000),  there  remain 
several  inherent  limitations  in  using  these  2D  cross-section-type  models.  Several  of  the  limitations 
associated  with  2D  modeling  can  be  eliminated  through  the  use  of  3D  models.  For  example,  2D 
models  cannot  account  for  variation  in  topography,  geology,  sources  or  sinks,  and  boundary 
conditions  along  the  length  of  the  soil  structure.  If  these  conditions  exist,  a  3D  model  may  be 
needed. 

Seepage-Based  Factor  of  Safety.  Beginning  with  a  laboratory  setup  as  shown  in  Figure  1, 
the  stability  of  the  soil  may  depend  on  whether  the  net  force  at  a  point  X  is  upward  or  downward. 
Suppose  the  soil  matrix  is  in  a  hydrostatic  condition,  then  the  upward  force  and  the  downward 
force  at  point  X  can  be  expressed  as 


Fx  = 


upward 


(A H  +  dponding  +L)jyw 


ponding 


(1) 


Fx 

rim 


downward 


(2) 


with  the  unit  weight  of  the  saturated  soil  ( ys )  defined  as 


y,  =  <I>Yw+(1-<I>)y^ 


(3) 
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where  AH  is  the  difference  in  head  across  the  soil  layer,  dponding  is  the  depth  of  water  above  the 
soil  matrix,  L  is  the  thickness  of  the  soil  matrix,  yw  is  the  unit  weight  of  water,  4>  is  the  porosity, 
and  ydry  is  the  unit  weight  of  the  dry  soil  matrix. 


Figure  1.  Factor  of  Safety 
determination. 


Three  scenarios  regarding  the  stability  of  the  soil  thus  exist: 

Critical  Condition:  When  F*pward  =  Fdownward ,  there  results 

(AH  +  dponding  +  L^jyw  =  dpondingyw  +  Lys 


which  can  be  rewritten  as 


AH 

~L~ 


(l-4>) 


Yd, 


ry 


=  (1-4>)(G,-1) 


where  Gs  is  the  specific  gravity  of  the  soil. 

Stable  Condition:  When  F*ward  <  Fdownward ,  there  results 


Unstable  Condition:  When  F*pward  >  F?ownward,  there  results 


AH 
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Therefore,  the  FoS  can  be  defined  as  the  ratio  of  the  stabilizing  downward  force  to  the  generally 
destabilizing  upward  force  applied  to  the  soil.  It  can  be  written  as 


FoS 


(1H>)(G,-1) 

(AH) 

y  L  , 


(8) 


where  icr  is  the  critical  hydraulic  gradient,  and  iv  is  the  exit  gradient  at  the  point  of  interest.  The 
critical,  stable,  and  unstable  conditions  are  thus  defined  as  FoS  =  1,  FoS  >  1,  and  FoS  <  1  (or 
icr  =  iv,  icr  >  iv,  and  icr  <  iv),  respectively.  Under  the  critical  condition,  the  vertical  effective  stress  is 
zero,  but  under  unstable  conditions,  the  potential  for  movement  of  soil  exists.  Due  to  variability  of 
factors,  the  critical  condition  is  unlikely  to  exist.  Therefore,  an  FoS  screening  criteria  >  1  is  chosen 
to  distinguish  stable  and  unstable  conditions  in  general  practice. 

In  general,  the  exit  gradient  at  a  point  of  interest  can  be  computed  from  the  simulated  head 
distribution  of  a  model  as  explained  in  Figure  2,  where  the  computational  model  accounts  for 
boundary  conditions,  sources/sinks,  any  overlying  material  at  the  point  of  interest,  and  the 
specified  ponding  depth  above  the  top  soil.  The  aforementioned  exit  gradient  across  the  soil  is 
assumed  to  be  constant  and  may  be  estimated  from  the  values  of  total  head  at  points  X  and  Y  and 
the  vertical  distance  between  the  points  (i.e.,  L  (Figure  2)).  The  total  heads  at  X  and  Y  are 
denoted  by  Hx  and  HY,  respectively,  in  Equation  9. 


AH   Hx-Hy 


(9) 


Figure  2.  Hydraulic  gradient  at  a  point-of-interest 
computation. 


According  to  the  U.S.  Army  Corps  of  Engineers  (USACE),  Jacksonville  District  Memorandum 
with  the  subject  Herbert  Hoover  Dike  Loadings  and  Factors  of  Safety  for  Use  in  Seepage  and 
Stability  Analyses  (U.S.  Army  Corps  of  Engineers  (USACE),  Jacksonville  District  2008),  the 
FoS  against  piping  of  foundation  soils  due  to  under  seepage  (FoSpipu,g)  and  the  FoS  against  uplift, 
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or  heave,  of  blanket  materials  (FoSupuft)  should  be  calculated  for  dike-safety  evaluation,  where 
FoSPiPing  is  for  pipeable  topsoil  (noncohesive  soils  such  as  sand)  and  FoSupHft  for  nonpipeable 
topsoil  (cohesive  soils  such  as  clay  and  silt).  The  FoS  for  pipeable  material  (i.e.,  FoSPiPing)  can  be 
computed  directly  using  Equation  8  while  FoSupiift  for  nonpipeable  materials  can  be  computed 
using  the  following  equation: 


FoS, 


2Xi  (Yz-Yw)^,- 


y  V"  A z. 


uplift 


AH 
A z 


A/Y 


(10) 


where  /?  is  the  number  of  soil  types  included  in  the  nonpipeable  layer,  c[).  is  the  porosity  of  the 
/-th  soil  type,  G,  is  the  specific  gravity  of  the  z-th  soil  type,  and  A z  —  Y-'f  A z{  denotes  the  total 
thickness  of  the  nonpipeable  layer. 

Note  that  Equations  8  and  10  used  to  compute  FoS  for  piping  and  uplift  are  essentially  the  same. 
However,  the  mechanism  for  initiation  is  slightly  different.  Piping  can  be  thought  of  as  a 
phenomenon  where  the  gradient  through  the  soil  is  such  that  movement  of  individual  soil  grains 
is  possible,  whereas  uplift  is  a  condition  when  the  pressure  beneath  the  blanket  exceeds  the 
effective  weight  of  the  blanket  material  and  water.  While  piping  deals  with  pipeable  materials  as 
sands,  uplift  considers  nonpipeable,  low-permeability  materials  such  as  clays  and  may  involve 
multiple  materials  of  different  thicknesses.  That  is,  the  computation  of  uplift  includes  the 
composite  effect  of  the  weight  of  all  the  layers,  but  the  piping  FoS  calculation  in  a  multilayer 
system  would  be  the  minimum  value  among  the  layers. 

METHODOLOGY  FOR  FACTOR  OF  SAFETY  COMPUTATION:  A  computer  program, 
written  in  Fortran  90,  was  developed  as  a  postprocessor  for  the  3D  WASH123D  groundwater 
model  (Yeh  et  al.  2006),  a  finite  element  computational  model  that  approximates  the  Richards 
equation  for  pressure  heads  at  the  nodes  of  unstructured  meshes.  The  postprocessing  utility  code 
computes  FoSPiping  using  Equation  8  with  the  exit  gradient  defined  as 


hb-ha 


AB 


(11) 


where  HA  is  the  total  head  at  a  boundary-face  node  A,  Hb  is  the  total  head  at  a  nonboundary- face 
node  B  that  is  associated  with  node  A,  and  dAB  is  the  distance  from  node  A  to  node  B. 

Generally,  the  distance  from  node  A  to  node  B  is  close  to  the  thickness  of  the  soil  associated 
with  node  A  divided  by  the  number  of  element  layers  used  to  discretize  the  material.  For 
instance,  if  the  soil  shown  in  Figure  3  associated  with  the  boundary-face  node  A  is  one  foot  (ft) 
thick  and  discretized  using  a  single  layer  of  elements,  then  the  distance  from  node  A  to  node  B 
should  be  approximately  a  foot. 
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Figure  3.  The  exit  gradient  computed 
at  boundary-face  node  A 
uses  the  total  heads  at  nodes 
A  and  B  and  the  distance 
from  nodes  A  to  B.  The  exit 
gradient  may  be  neither 
vertical  nor  horizontal. 

The  utility  code  requires  the  following  data  as  input: 

•  The  mesh  of  the  3D  groundwater  model  of  interest,  which  includes 

o  nodal  coordinates  and 
o  element  indices. 

•  The  head  distribution  over  the  entire  domain  of  the  3D  model  (total  or  pressure  head). 

•  The  parameters  necessary  to  compute  FoS,  which  includes 

o  specific  weight  of  water, 
o  specific  weight  of  soil  materials  (solids), 
o  porosity  of  soil  materials,  and 
o  FoS  screening  criteria. 

The  utility  code  computes  the  critical  gradient  and  the  FoS  value  associated  with  each  of  the 
boundary- face  nodes  of  interest.  To  help  the  analysis,  it  also  outputs  the  exit  gradient  computed 
at  each  node  of  interest. 

Computation  of  Exit  Gradient  at  Unstructured  Mesh  Nodes.  In  the  preceding,  the  FoS 
and  the  exit  gradient  were  considered  only  in  the  vertical  direction.  O’Leary  et  al.  (2013) 
provided  an  overview  of  the  available  methods  that  can  be  used  to  evaluate  critical  horizontal 
seepage  gradients.  The  angle  of  internal  friction  that  can  be  determined  in  the  laboratory  by  a 
direct  shear  test  or  a  triaxial  stress  test  needs  to  be  taken  into  account  when  the  horizontal  FoS  is 
evaluated.  Al-Hussaini  et  al.  (1997)  proposed  to  use  the  following  equation  to  compute  critical 
horizontal  gradient: 
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• horizontal  •  /«1  o\ 

=  icr  tan  cp  (12) 

where  cp  is  the  angle  of  internal  friction.  The  angle  of  internal  friction  ranges  from  20°  (clay)  to 
48°  (gravel  with  some  sand).  As  a  result,  the  coefficient  of  internal  friction  (tancp )  ranges  from 
0.364  to  1.11,  making  the  critical  gradient  either  lower  or  higher  than  the  vertical  gradient.  While 
the  utility  uses  Equation  12  to  compute  the  horizontal  critical  gradient,  other  methods  can  also  be 
incorporated  into  the  utility  for  better  estimates. 

Due  either  to  the  topographical  variation  or  to  the  discretization  of  unstructured  mesh,  the  exit 
gradient  computed  can  be  in  a  direction  that  is  neither  vertical  nor  horizontal  (Figure  3).  To 
handle  this  situation,  the  utility  code  computes  both  vertical  FoS  and  horizontal  FoS  using  the 
computed  head  at  the  nonboundary- face  node  closest  to  the  boundary- face  node  of  interest.  For 
example,  node  B  is  the  nonboundary- face  node  closest  to  the  boundary-face  node  A  in  Figure  3. 
The  vertical  FoS  (i.e.,  FoSA™cal )  and  the  horizontal  FoS  (i.e.,  FoShABzontal)  associated  with  nodes 
A  and  B  can  thus  be  computed  as 


FoS 


o vertical  ^ cr  _ 

i 

bo 

r°^AB  ~  ~  1 

l 

\ 

hb-ha  | 

d AB 

"  (h,-ha) 

—  C"  (G,-l)tan<!  ( G,-l)dAB 

AB 

l 

[hb-ha 

K  d AB 

(H,-Ha) 

tan  <p 


(13) 


(14) 


where  HA  is  the  total  head  at  a  boundary-face  node  A,  HB  is  the  total  head  at  a  nonboundary-face 
node  B,  and  cIAb  is  the  distance  from  node  A  to  node  B. 

With  the  fact  that  only  the  vertical  FoS  needs  to  be  evaluated  when  the  boundary  face  is 
horizontal,  and  only  the  horizontal  FoS  is  needed  when  the  boundary  face  is  vertical,  the  final 
FoS  for  a  location  on  a  sloped  boundary  face  can  be  estimated  using  a  simple  equation 
accounting  for  the  contributions  from  both  the  vertical  FoS  and  the  horizontal  FoS  as 

FoSpiping  =  FoSZtiml  cos  a  +  FoShA°Brizontal  sin  a  (15) 


where  cos  a= — — ,  sina  =  — — ,  and  A zAB  and  Axab  are  the  vertical  and  the  horizontal 

d ab  dAs 

distances,  respectively,  between  nodes  A  and  B. 

To  calculate  FoS  adequately,  the  nonboundary-face  node  used  should  be  placed  in  a  location 
such  that  the  straight  line  connecting  it  and  the  respective  boundary-face  node  is  approximately 
perpendicular  to  the  sloped  boundary  face.  This  needs  to  be  taken  into  account  during  the 
development  of  the  computational  mesh  of  the  3D  groundwater  model. 
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While  the  utility  provides  a  quick  estimate  of  FoSpiping  at  each  of  the  boundary-face  nodes  of 
interest,  more  field  data  must  be  incorporated  to  compute  FoSupuft  when  non-pipeable  surface 
materials  (e.g.,  clay  blanket)  are  present.  Moreover,  although  this  utility  provides  an  effective 
way  to  identify  the  potential  for  initiation  of  internal  erosion,  additional  considerations  should  be 
given  to  the  likelihood  that  backward  erosion  piping  would  continue  towards  the  reservoir  and 
cause  a  breach.  Additional  soil  mechanics  information  must  be  taken  into  account  when 
completing  internal  erosion  studies. 

DEMONSTRATION  EXAMPLE:  A  WASH123D  model  was  developed  within  Reach  2  of  the 
Herbert  Hoover  Dike  (HHD)  around  Lake  Okeechobee,  Florida,  to  study  the  feasibility  of  a 
proposed  Seepage  Management  Test  Facility  (SMTF)  and  a  shallow  cutoff  wall  alternative  to 
mitigate  seepage,  as  detailed  in  the  HHD  Flow  Model  Report  (Pickett  et  al.  2013).  Figure  4 
depicts  the  typical  embankment  of  the  SMTF  taken  from  construction  plans.  The  HHD  was  built 
to  prevent  flooding  of  Lake  Okeechobee  and  allowed  for  the  development  of  the  Everglades 
Agricultural  Area  (EAA).  The  SMTF  includes  a  drainage  system  built  into  the  downstream 
slope  of  the  HHD  and  includes  an  inclined  chimney  drain  to  collect  through-embankment 
seepage,  a  trench  drain  to  intercept  through-foundation  seepage  for  pressure  relief  and  a 
blanket  drain  to  convey  the  seepage  collected  by  the  chimney  and  trench  drains  to  the  ditch. 
The  shallow  cutoff  wall  consists  of  an  impermeable  material  installed  along  the  HHD  crest  to 
inhibit  horizontal  seepage  through  the  dike. 


LAKE  OKEECHOBEE 


Figure  4.  Typical  embankment  and  instrumentation  cross  section  of  the  SMTF  from  construction  plans 
provided  by  the  USACE,  Jacksonville  District. 


Model  Discretization.  The  model  domain  covers  an  area  of  0.29  square  miles  and  was 
discretized  with  163,884  nodes  and  570,940  elements.  The  mesh  resolution  varies  to  represent  all 
of  the  SMTF  and  relevant  surrounding  features.  This  model  consists  of  multiple  element  types 
with  tetrahedrons  on  top  of  triangular  prisms  (Figure  5)  to  allow  pinch  out  of  material  regions 
and  additional  resolution  of  the  SMTF  drain  in  the  Z-axis  direction. 


7 


ERDC/CHL  CHETN-XI-2 
August  2014 


Materials 


Undifferentiated  Fill 
Peat  {  Silt 
Interbedded  Zone 
Zone  1 

Confining  Unit  1 
Zone  2 

Confining  Unit  2 
Zone  2 
Sand 

Cutoff  Wall 
Fine  Drainage  Gravel 
Coarse  Drainage  Gravel 
New  Embank ment  Fill 


View 

Lx 


Figure  5.  The  WASH123D-SMTF  model  mesh  with  material  assignments  (full  domain,  top;  close 
up  of  dike,  bottom)  shown  vertically  stretched  (z-magnification  =  5). 

Geology  and  Configuration.  As  shown  in  Figure  6,  the  drainage  system  stretches  three- 
quarters  of  the  width  of  the  model  domain  and  two-thirds  of  the  landside  HHD  face  (from  toe 
ditch  to  crest)  with  the  upper  chimney  resting  in  the  Undifferentiated  Fill,  covered  by  New 
Embankment  Fill,  and  the  lower  chimney  cutting  through  the  Peat/Silt  and  Interbedded  Zone 
regions  into  Zone  1 .  The  New  Embankment  Fill  overlays  a  larger  area  than  the  Drainage  Gravels 
due  to  construction  grading  (Figure  4)  and  is  assumed  to  retain  original  Undifferentiated  Fill 
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material  properties  from  proper  compaction  (materials  are  distinguishable  only  by  name).  The 
drain  face  is  approximately  3  ft  above  the  toe  ditch  thalweg  that  incises  into  the  Peat/Silt  region. 
Figure  6  also  shows  the  shallow  cutoff  wall  configuration,  which  extends  the  entire  width  of  the 
domain  and  terminates  at  the  interface  of  the  Interbedded  Zone  and  Zone  1.  Figure  7  presents  the 
four  configurations  considered  in  the  study:  existing,  shallow  cutoff  wall,  SMTF  drainage  system 
(hereafter  referred  to  as  just  drain),  and  shallow  cutoff  wall  with  drain. 

Simulation  Results.  Five  sets  of  boundary  conditions  were  considered  in  the  study.  However, 
only  the  set  with  an  extremely  high  lake  stage  (i.e.,  30  ft  NAVD88)  is  presented  herein  for 
demonstration.  Figure  8  shows  the  computed  total  head  distributions  associated  with  the  four 
configurations  from  WASH123D.  Figure  9  depicts  the  areal  extents  where  the  initiation  of  piping 
would  likely  occur  due  to  high  hydraulic  gradients  according  to  the  FoS  calculated  using  the 
utility  code.  The  FoS  analysis  was  conducted  on  the  entire  surface  nodes  of  the  model.  To  filter 
the  results,  four  values  (1,  1.25,  1.5,  and  2)  were  used  as  screening  criteria.  It  was  assumed  that 
areas  with  FoS  values  >  2  would  always  be  stable. 

Focusing  on  the  Peat/Silt  material  region  of  the  toe  ditch,  the  shallow  cutoff  wall  provides  nearly 
no  protection  because  it  does  not  regulate  flow  into  the  toe  ditch.  It  is  effective,  though,  in 
increasing  the  FoS  on  the  embankment  surface  (not  shown  in  Figure  9  since  FoS  >  2). 
Conversely,  the  drain  effectively  reduces  the  hydraulic  gradient  at  the  toe  ditch  bottom  by 
diverting  the  flow  and  increasing  the  FoS  as  a  result.  The  combination  of  the  shallow  cutoff  wall 
and  drain  had  negligible  gain  from  the  drain  alone.  This  result  is  better  seen  in  Tables  1  and  2. 


Materials 


Undifferentiated  Fill 
Peat  /  Silt 
Interbedded  Zone 
Zone  1 

Confining  Unit  1 
Zone  2 

Confining  Unit  2 
Zone  3 
Sand 
Cutoff  Wall 
Fine  Drainage  Gravel 
Coarse  Drainage  Gravel 
New  Embankment  Fill 


View 

u 


- Planar  Drain  Extents 


Figure  6.  The  WASH123D-SMTF  model  geology,  shallow  cutoff  wall,  and  drain  regions. 
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Figure  7.  The  WASH123D-SMTF  model  configurations. 
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Figure  8.  The  WASH123D-SMTF  model  total  head  distribution  along  the  embankment. 
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Figure  9.  The  WASH123D-SMTF  model  FoS  distribution  along  the  embankment  (close  up  of 
HHD  toe  ditch). 


Table  1.  FoS  Distribution  on  WASH123D-SMTF  Model  Surface  Nodes 

Value 

Number  of  Occurrences  i 

'per  model  configuration) 

Existing 

Shallow  Cutoff 
Wall 

Drain 

Shallow  Cutoff 
Wall  &  Drain 

FoS  <  1.00 

681 

681 

423 

422 

1.00  <  FoS  <  1.25 

7 

7 

20 

20 

1 .25  <  FoS  <  1 . 50 

0 

0 

47 

43 

1.50  <  FoS  <2.00 

2 

2 

101 

99 

2.00  <  FoS 

1,958 

1,810 

1,167 

1,120 

Unquantifiable 

3,902 

4,050 

4,792 

4,846 

Note:  FoS  values  were  calculated  at  all  quantifiable  surface  locations  where  hydraulic  gradient 
caused  water  to  exit  the  ground.  Surface  mesh  of  6,550  total  nodes,  HHD  toe  ditch  thalweg 
consists  of  201  nodes  with  an  approximately  10  ft  spacing,  and  EAA  minor  ditch  consists  of  68 
nodes  with  an  approximately  30  ft  spacing. 
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Table  2.  Comparsion  of  Average  FoS  of  WASH123D-SMTF  Model  Configurations 

Configuration 

Average  FoS 

Average  FoS  Increase  %  (from  Existing) 

Existing 

0.215 

Shallow  Cutoff  Wall 

0.219 

1 .7% 

Drain 

0.755 

250.7% 

Shallow  Cutoff  Wall  &  Drain 

0.749 

247.7% 

Note:  Only  surface  locations  with  FoS  between  0  and  2  were  used  to  calculate  average  values. 
Change  in  percentage  values  are  the  increase  in  FoS  caused  by  the  proposed  configuration 
compared  to  the  existing  configuration 


Table  1  lists  the  number  of  surface  locations  within  various  FoS  criteria  ranges.  Table  2  lists  an 
overall  average  FoS,  when  considering  just  surface  locations  with  FoS  <  2,  and  percent  change  in 
the  overall  average  FoS  comparing  the  configurations  to  the  respective  existing  condition  of  each 
boundary  condition  set.  These  occurring  failure  results  rely  on  a  regional  geology  conceptual 
model  that  simplifies  heterogeneity  local  to  the  HHD  and  underestimates  existing  conditions 
preventing  failure,  such  as  vegetation.  However,  note  that  the  change  of  FoS  due  to  controlling 
seepage  is  a  legitimate  and  usable  result. 

SUMMARY:  A  utility  code  was  developed  to  compute  FoSPiping  at  boundary-face  nodes  using 
the  computed  heads  from  3D  groundwater  simulations.  The  utility  code  computes  exit  gradients 
and  uses  a  simple  equation  (Equation  14)  to  estimate  the  FoSPiping  at  locations  on  a  sloped 
boundary  face  based  on  both  the  vertical  FoS  (Equation  12)  and  the  horizontal  FoS  (Equation 
13).  A  groundwater  model  developed  for  an  SMTF  feasibility  study  was  employed  for 
demonstration.  The  model  accounts  for  a  drainage  system  as  well  as  a  cutoff  wall.  Four 
configurations  were  compared  to  help  identify  effective  components  for  seepage  management. 
Note  that  the  FoS  results  of  the  model  rely  on  a  regional  geology  conceptual  model  that 
simplifies  heterogeneity  local  to  the  HHD.  Future  evaluations  of  the  WASH123D-SMTF  model 
would  benefit  from  the  incorporation  of  epistemic  uncertainty.  Additional  knowledge  and  data 
should  necessitate  increasing  the  mesh  resolution  that  would  refine  the  FoS  analysis.  For  future 
improvement,  Equation  11  can  be  incorporated  into  the  utility  to  compute  FoSupuft  for 
nonpipeable  surface  materials. 

ADDITIONAL  INFORMATION:  This  work  was  completed  for  an  HHD  groundwater  modeling 
project,  funded  by  the  U.S.  Army  Corps  of  Engineers  (USACE),  Jacksonville  District.  Dr.  Stacy 
E.  Howington  of  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Coastal 
and  Hydraulics  Laboratory  (CHL);  Laura  D.  Bittner,  P.E.,  Chief  of  Hydrology,  Hydraulics,  and 
Coastal  Branch,  of  the  USACE,  Philadelphia  District;  and  Bradley  Cox  of  the  USACE, 
Jacksonville  District,  reviewed  this  technical  note  and  provided  valuable  comments.  For 
additional  information,  contact  Hwai-Ping  (Pearce)  Cheng,  ERDC/CHL,  3909  Halls  Ferry  Road, 
Vicksburg,  MS  39180,  at  601-634-3699  or  e-mail:  Hwai-Pina .Chenci (a)usace.  grrnujnil ■  This 
CHETN  should  be  cited  as  follows: 

Cheng,  H.-P.,  K.  D.  Winters,  S.  M.  England,  and  R.  E.  Pickett.  2014.  Factor  of 

safety  analysis  using  3D  groundwater  simulation  results.  ERDC/CHL  CHETN- 

XI-2.  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development  Center. 
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An  electronic  copy  of  this  CHETN  is  available  from  /iffp://c/t/.er(ic.usace.army.mi l/chetn. 
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